UNCLASSIFIED 


AD  633  718 

THE  GENERALIZED  CA YL E Y -HAM  I L TON  THEOREM 
IN  N  DIMENSIONS 

John  S.  Lew 

Brown  Un i vers i ty 
Providence,  Rhode  Island 

March  4,  1965 
Processed  for .  .  . 

DEFENSE  DOCUMENTATION  CENTER 
DEFENSE  SUPPLY  AGENCY 


U.  *  DEPARTMENT  OF  COMMERCE  /  NATIONAL  BUREAU  OF  STANDARDS  /  INSTITUTE  FOR  APPLIED  TECHNOLOGY 


UNCLASSIFIED 


NASA  Scientific  and  Technical  Information  Facility 

operated  for  the  National  Aeronautics  and  Soace  Admimst’ahon  by  Documentation  Incorporated 


/ 


Post  Office  Box  33 
College  Park  Md  20740 


Telephone  Area  Code  301 
|  7792121 

.  X- 

FACILITY ’CONTROL  NO.  _  ..  C 


DATE  ' 


•  '><  x 
^  &  C 


ATTACHED  IS  A  DOCUHENT  ON  LOAN 


NASA  Scientific  and  Technical  Information  Facility 


TO:  Defense  Documentation  Center 

Attn:  DDC-IRC  (Control  Branch) 
rameron  Station 
Alexandria,  Va.  223IA 


in  accordance  with  he  NASA-DOD  Cooperative  AD  Number  Assignment  Agreement  it  is 
requested  that  an  AD  number  be  assigned  to  the  attached  report. 


Br, 


As  this  is  our  only  available  copy  the  return  of  the  document  (with  AD 
number  and  any  applicable  distribution  I i notations)  to  the  address  below 
is  essential . 


□  This  document  may  be  retained  by  DDC.  If  retained,  please  indicate  AD 
number  'and  any  applicable  distribution  limitations  on  the  reproduced 
copy  of  the  tit'e  page  and  return  to  the  address  below. 


Return  Address:  NASA  Scientific  and  Technical  Information  Facility 
Attention:  INPUT  BRANCH 
P.  0.  Box  33 

College  Park,  Maryland  207^0 


FFNo  2M*  Revised  Aug  65 


AD  633? 1 g 


Division 
►  of 


Applied  Mathematics 


ONR  Technical  Repeat  No»  6 
ARPA  Technical  Report  No*  AM»20 

(^THE  GENERALIZE?  CAYLEY-HAMILTON  THEOREM 
IN  n  DIMENSIONS 


by 

John  S*  Lew 


Brown  University 


Providence 
Rhode  Island 


The  Generalized  Cayley-Hamilton  Theorem  in  n  Dimensions 

John  S.  Lew 
Brown  University 

In  1945,  Reiner  ,  by  means  of  the  Cayley-Hamilton  theorem, 
obtained  a  canonical  form  for  a  polynomial  relation  between  a 
stress  matrix  and  a  strain-velocity  matrix^  since  that  time  an 
extensive  theory  has  been  developed  for  canonical  forms  of  non¬ 
linear  constitutive  equations.  More  recently,  for  a  polynomial 
relation  between  ne  tensor  and  a  number  of  other  tensors,  the 
problem  of  finding  the  restrictions  imposed  by  a  symmetry  group 
was  reduced  by  Smith  and  Riviin^,  and  by  Pipkin  and  Rivlin^ , 
to  that  of  finding  an  integrity  basis  for  a  set  of  such  tensors ; 

and  then,  for  the  full  (or  proper)  orthogonal  group  in  Euclidean 

(4 ) 

2-space  or  3-space,  such  a  basis  was  determined  by  Rivlinv  , 

(tt)  (&) 

Spencer  and  Rivlinv  ,  and  Spencer  ,  and  its  irreducibility 

( 7} 

proven  by  Smith  • 

In  this  development,  an  important  tool  has  been  a  generaliza¬ 
tion  of  the  Cayley-Hamilton  theorem,  in  2-space  or  3-space,  from 
one  to  several  matrix  variables During  this  time,  it  has 
been  clear  that  the  corresponding  identity  in  n-space,  for  any 
particular  n,  could  be  obtained  i'**  a  finite  but  discouraging 
number  of  steps j  however  the  form  of  this  relation  for  an  arbi¬ 
trary  n  has  not  been  given.  We  shall  obtain  this  form,  which  is 
the  intuitive  generalization  of  the  results  in  two  and  three 

(1)  See  Reference  4-. 

(2)  See  Reference  7. 

(3)  See  Reference  3» 

(4)  See  Reference  5« 

(5)  See  References  9,10,11. 

(6)  See  Reference  8* 

(7)  See  Reference  6. 

(8)  See  Reference  5« 
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dimensions,  and  note  that,  properly  expressed,  it  is  a  polynomial 
relation  in  the  given  matrices  and  their  traces  with  all  coeffici¬ 
ents  +1. 

For  an  arbitrary  real  or  complex  nxn  matrix  A  the  Cayley- 
Hamilton  theorem  states  that 

1)  s  (-dV^s.Ca)  =  0 

i=o  1 

where  A0  =  I,  and  s^A)  is  the  i*th  symmetric  polynomial  in  the 
characteristic  roots  of  A*  If  we  let  t^CA)  =  tr  A^  for 

(Q) 

j  =  1,2, •• •  then  the  well-known  relations 

S1  =  tl 

2)  2s 2  =  s^  -  t2 

3s3  =  s2t1  -  sxt2  +  t^ 

and  so  forth  can  be  solved  recursively  for  each  si  in  terms  of 
t1? •• • ,  ti  to  yield 

lls1  =  ^ 

3)  2!s2  =  -  t2 

3's,  =  *1  "  ’tit2  +  2t3 

and  so  forth.  Thus  the  Cayley-Hamilton  theorem  can  be  expressed 
as  a  relation  in  A  and  the  t j (A) • 

In  one  dimension  this  process  gives 

k)  A  -  I  tr  A  =  0 

which  is  trivial;  and  in  two  dimensions  it  gives 

5)  A2  -  A  tr  A  +  |l[(tr  A)2  -  tr  A2]  =  0. 


(9)  See  p.  9  of  Reference  12. 


that 


If  we  apply  to  this  equation  the  polarization  operator  dg^, 

is,  if  we  replace  A  by  A+xB,  for  a  real  variable  x,  and  evaluate 

the  derivative  in  x  at  the  point  x  =  0,  then  we  obtain 

6)  AB  +  BA  -  A  tr  B  -  B  tr  A  +  l[tr  A  tr  B  -  tr  AB]  =  0 

which  is  the  generalized  identity  in  two  dimensions'0^#  Note 
in  this  result  that  all  permutations  of  A  and  B  appear,  since 
A  and  B  need  not  commute,  but  that  the  fraction  ^  disappears, 
since  scalars  commute  and  traces  have  cyclical  symmetry. 

In  three  dimensions  equations  (1)  and  (3)  give 

7)  A3-  A2tr  A  +  |Af(tr  A)2  -  tr  A2] 

-  |l[(tr  A)3  -  3  tr  A  tr  A2  +  2  tr  A3]  =  0 

applied  to  which  the  polarization  operator  dg^  again  yields  a 
relation  in  two  variables.  However,  we  desire  a  completely 
polarized  relation,  in  which  no  matrix  variable  has  degree  more 
than  unity.  Thus  if  we  also  apply  d^  for  another  3x3  matrix  C, 
and  lets  denote  the  sum  over  all  permutations  of  (A,B,C),  then 
we  obtain 

8)  0  =  SABC  -  SAB  tr  C  +  SA[tr  B  tr  C  -  tr  BC] 

-  I[tr  A  tr  B  tr  C  -  tr  A  tr  BC  -  tr  B  tr  CA  -  tr  C  tr  AB 
+  tr  ABC  +  tr  CBA] 

(11) 

which  is  the  generalized  identity  in  three  dimensions'  '•  Note 
here  again  that  all  fractions  disappear  by  the  properties  of 
scalars  and  traces. 


(10)  See  Reference  5« 

(11)  See  Reference  5* 


Now  in  the  derived  expression  for  each  s^ ,  the  terms  corres 
pend:  to  partitions  of  i,  that  is,'  to' sequences .p.  =  of 
non-negative  integers  with  q(p)  =  i,  where 

CO  CO 


9)  p(p)  =  2  (j-l)m.  ,  q(|i)  =  2  jm,  . 

j=l  3  j=l  3 

Each  is  interpreted  as  the  number  of  subsets  containing  pre- 

v 

cisely  j  elements  in  the  corresponding  subdivision  of  a  set  con¬ 
taining  precisely  i  elements,  so  that  clearly  *  =  0  for  j  >  i  and 

J 

thus  such  sequences  have  all  entries  but  a  finite  number  equal  to 
zero.  If  we  let  -s  denote  the  sequence  (t-j,  ,t0,,..)  with  t ^  =  tr 
as  before,  then  we  may  let 


m-,  nu 

10)  ...  , 

a  product  which  thus  has  all  factors  but  a  finite  number  equal  to 
unity. 


Each  partition  p  with  q(p)  "  q  labels  a  conjugate  class  C 
in  the  group  S  containing  all  permutations  of  q  elements,  namely 

H 

that  class  in  which  all  permutations  may  be  factored  into  disjoint 

cycles  of  which  have  length  1,  m2  have  length  2,  and  so  forth. 

The  parity  of  all  elements  in  Cn  is  easily  shown  to  be  sgn(p)  = 

r 

(_1)P^) ,  and  the  number  of  elements  in  is  well-known  to  be^12^ 

m,  mp 

11)  c(p)  =  ql/m^Jm,,!...!  x2  ... 


a  quotient  whose  denominator  has  all  factors  but  a  finite  number 

equal  to  unity.  However,  the  general  expression  of  the  set  (3) 

(1^) 

can  then  be  written v 


12)  iJs.  =  Z  sgn(p)c(p>rP 

1  q(p)=i 


(12)  See  3.6  of  Reference  1  or  IV. 4  of  Reference  2. 

(13)  See  6.2  of  Reference  1  or  (4.24)  of  Reference  2. 


and  the  form  (1)  of  the  Cayley-Hamilton  theorem  can  be  rewritten 

13)  s  (-l)V'1  Z  sgn(!.)c(ii)T,l(A)/iJ  =  C. 

i=o  q  (p,  )~i 

Now  we  need  only  completely  polarize  this  equation,  noting 
that  it  is  homogeneous  of  degree  n  in  A;  that  is,  we  need  only 
replace  the  n  equal  variables  A  by  all  permutations  of  n  distinct 
variables  A^,*.«,An,  and  put  the  sum  of  all  such  expressions 
equal  to  zero.  For  each  i  the  corresponding  term  in  (13)  then 
yields  ni  terms,  of  which  we  may  collect  all  those  terms  such 
that  \(i)  ?  •  ••  )  >  In  any  order,  appear  in  the  inner  sum,  and 

Atc(1+1)  ’  ••  *  ,Aft(n)  ’  in  any  order?  appear  in  the  outer  sum.  Thus 
for  each  i  we  obtain  a  sum  over  the  (^)  ways  to  select  a  subset 
of  i  elements  from  {A^,,.,  ,Ar} ,  with  each  summand  of  the  form 

14) 

( —1 )  [2  all  permutations  of  A^  ^  *  *  *A7t(u)  aef f .  (A^^j^  >  ••  •  > i )  ^ 

and  to  find  the  coefficient  for  each  selection  we  need  only 
replace  the  i  equal  variables  A  by  all  permutations  of  i  distinct 
variables  A„ (]_)>••• 

m,  j  Di  m,  fflj 

15)  2  sgn(p.)(tr  A)  A,..(tr  A1)  /m, l...m,  11  ...i 

q (p )=i  1  1 

and  take  the  sum  of  all  such  expressions. 

But  for  each  p,  in  the  sum  (15)  many  of  the  resulting  ii  terms 

are  equal 5  in  particular,  for  each  3  the  m^  traces  of  products 

containing  j  factors  may  be  permuted  in  all  m^J  ways  without 

changing  the  result,  and  in  each  of  these  traces  the  j  factors 

may  be  cycled  in  j  ways  without  changing  the  result.  Since  the 

various  terms  are  otherwise  distinct,  the  order  of  their 


-1 
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degeneracy  is  precisely  the  denominator  associated  with  the 
given ji  in  the  sum  (15),  and  thus  the  coefficient  is  precisely 
the  siim,  over  all  partitions  of  i  and  all  essentially  distinct 
permutations  of  the  A^.^,  of  terms 
mi  m2 

16)  BgnOO nitp(A,|(Bi+2j_1)lAK(Bi+j,j))...  . 

In  summary,  the  generalized  Cayley-Hamilton  theorem  in  n 
dimensions  asserts  the  vanishing  of  the  sum  for  i  =  0,...,n  of 
a'1  L  essentially  distinct  terms  of  the  form  (l1*),  in  which  the 
coefficient  is  the  sum  for  q(p.)  =  i  of  all  essentially  distinct 
terms  of  the  form  (16).  Furthermore,  the  coefficients  are 
independent  of  n,  and  have  the  forms  given  in  equation  (8)  for 
I  =  1>2,35  finally,  oy  the  principle  just  stated,  the 
coefficient  for  i  =  k  has  the  form 

17) 

trAtrBtrCtrD  -  trAtrBtrCD  -  trAtrCtrBD  -  trAtrDtrBC  -  trBtrCtrAD 
-  trBtrDtrAC  -  trCtrDtrAB  +  trA(trBCD+  trDCB)  +  trB(trACD+  trDCA) 
+  trC  (trABD  +  trDBA)  +  trD(trABC+  trCBA)  +*.  trAB-irCD  +  trACtrBD 
+  trADtrBC  -  trABCD  -  trABD C  -  trACBD  -  trACDB  -  trADBC  -  trADCB 

Since  the  polarization  process  d^  can  be  defined  over  any 

( 1 M 

field  of  characteristic  zero  ,  these  results  are  all  valid 
over  any  such  field.  Indeed  since  the  coefficients  are  simpler 
in  the  completely  polarized  equations,  these  results  xnay  well  be 
provable  directly  In  n  variables,  rather  than  through  (13)* 
However,  this  discussion  indicates  the  explicit  form  of  the 
desired  relation,  and  thus  reduces  the  labor  of  deriving  it  to 
merely  that  of  writing  it  down. 


(I1*)  See  p.  of  Reference  12. 
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